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Abstract 

We investigate the phase transitions characterized by deconfinement and restoration of chiral and 
axial symmetries, at finite temperature, in the framework of QCD inspired models. We compare 
the results obtained in the SU(2) Polyakov-Nambu-Jona-Lasinio (PNJL) model with anomaly and 
" in its extended version, the Entangled Polyakov-Nambu-Jona-Lasinio (EPNJL) model. In the last 

q; version four-quark vertices with entanglement between the chiral condensate and the Polyakov loop 

are considered. The thermodynamics of the phase transitions, the meson spectrum, in particular the 
convergence of axial and chiral partners will be analyzed, as well as the topological susceptibility. 
We find that an explicit temperature dependence of the coupling vertices is necessary in both 
models in order to have effective restoration of the Ua(1) symmetry. 
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I. INTRODUCTION 



The study of matter under extreme conditions of temperatures and/or densities is cur- 
rently an important field of research in physics both experimentally and theoretically. In 
this limit, quantum chromodynamics (QCD), the theory of strong interactions, predicts that 
matter becomes a plasma of deconfmed quarks and gluons. In spite of the success of many 
aspects of the strong interaction physics, there are important features to clarify, such as 
the behavior of the masses of the mesons and the thermodynamics of strongly interacting 
matter in the region of the transition from hadronic to quark matter. 

Monte-Carlo simulations in lattice gauge theory provide a powerful non-perturbative first 
principle approach to QCD, but the region of large densities and low temperatures essen- 
tially remains inaccessible to lattice simulations. This encourages the study and enlargement 
of effective models with QCD symmetries such as Polyakov-Nambu-Jona-Lasinio (PNJL) 
model. The PNJL model is an effective model which respects important symmetries of QCD 
action. It contains quarks as fundamental degrees of freedom allowing for a self-consistent 
description of chiral symmetry breaking, a key feature of QCD in its low temperature and 
density regime. Besides, the coupling to the Polyakov loop allows to describe the (statis- 
tical) confinement/deconfinement phase transition by taking into account a static gluonic 
field in which quarks propagate jl-4]. The NJL part of the model (where a point-like inter- 
action mimicking frozen gluon is introduced) can describe the chiral phase transition and 
its associate order parameter, the quark condensate. On the other hand, it is known that 
in pure gauge theory there is a confined/deconfined phase transition whose associate pa- 
rameter is the Polyakov loop. Quarks couple simultaneously to the chiral condensate and 
to the Polyakov loop, allowing to examine the correlation between confinement and chiral 
symmetry breaking which direct relation is not yet shown in QCD. 

The PNJL model has proven to be successful in reproducing lattice data concerning QCD 
thermodynamics jsj]: the coupling to the Polyakov loop, resulting in a statistical suppression 
of the unwanted quark contributions to the thermodynamics below the critical temperature, 
allows a good description compatible with lattice results. It has been also used to describe 
the full phase diagram, essentially at the mean field level. 

An important query on QCD thermodynamics is the proximity or coincidence of the two 
phase transitions characterized respectively by restoration of chiral symmetry and decon- 



finement. Lattice QCD results concerning this subject have been a matter of debate and, 
while for Nt = 2+1 flavors these two phase transitions are reported to take place at dis- 

n 

tinct temperatures I5J, for Nt = 2 flavors results indicate that the two phase transitions 

n 

occur at the same temperature |6[, Tq = 174(3) (6) MeV. This effect could be the result of 
strong correlations (entanglement) between the quark condensate and the Polyakov field, <3>. 
An extension of the PNJL model has been proposed where this entanglement is taken into 
account, by endowing the four quark interaction of the model with a dependence on the 
°olyakov field. This is the so-called Entangled Polyakov-Nambu-Jona-Lasinio (EPNJL) 
3, In this model the two phase transitions occur at the same temperature and it is 
possible to reproduce lattice values. 

A longstanding question is whether the Ua(1) symmetry, that is broken in the vacuum at 
the quantum level by instantons, is still broken in the chiral symmetric phase. If the amount 
of the Ua(1) symmetry breaking decreases with temperature, the question is what fraction of 
it remains above the critical temperature and whether and when this symmetry is restored. 
It has been pointed out that there are phenomenological consequences for the nature of the 
ohase transition depending on the degree of anomaly present at the critical temperature 



10|. Several observables can exhibit signals of the restoration of axial symmetry, like the 



topological susceptibility, the meson axial chiral partners and the rf mass. Lattice calcula- 



tions have found evidence of the decrease of the topo 
in former works 



qgical susceptibility with temperature 



ll| as well as in recent ones [9|, llfj, [12 ] . In the last works attention has 
also been given to the meson correlators of chiral and axial partners, which should become 
degenerate when both symmetries are restored. Concernin g the rf , the decrease of its mass 



in medium has been predicted in several theoretical works |13Ml6|. and recent experimental 
results [l?! in Au+Au collisions are compatible with a decrease of about 200 MeV, which 
could indicate the return of the 9th "prodigal" Goldstone boson. 

QCD inspired models have also been used to study the restoration of axial symmetry, in 



particular the NJL model 



15 



18] and its extended version, the PNJL model 19j, |20|. A 



procedure used in several works to account for the decrease of the top qlog ica 
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susceptibility 
. In previous 



is to allow a temperature dependence of the anomaly coefficient 
works, we have discussed the possible relation between the restorations of chiral and axial 
symmetries, both in SU(3) NJL model and PNJL models including the 't Hooft interaction 



2l| . It has been found that, by allowing the presence of high momentum quarks at finite 



temperature (using an infinite cutoff as explained in 2l|) observables related with the Ua(1) 
symmetry breaking vanish (the topological susceptibility as well as the mass difference be- 
tween the meson axial partners) as a natural consequence of the effective restoration of chiral 
symmetry, without the need of additional assumptions. However, the situation is different 
in the framework of the SU(2) NJL and PNJL models, as we have shown in a recent work 



22|. In this case, an additional mechanism that insures the suppression of instantons with 



temperature is needed in order that the meson axial partners become degenerate. This was 
achieved by taking the anomaly coefficient as a decreasing function of the temperature. If 
such a mechanism is not considered, and although the meson chiral partners degenerate 
and the topological susceptibility vanishes, signals of the anomaly will remain long after the 
critical temperature since the meson axial partners do not converge. 

The goal of the present paper, besides enlarging the previous study exploring different 
schemes of temperature dependence of the coupling coefficients, is to extend it to the SU(2) 
EPNJL model. A comparative study will be performed and the question whether there is 
any relation between the deconfinement, restoration of chiral and axial symmetries will be 
analyzed. Since restoration of axial symmetry in the PNJL model can only be achieved 
by assuming a temperature dependence of coupling coefficients, we will discuss whether an 
explicit temperature dependence of the vertices is also necessary in the EPNJL model for 
this purpose or if the temperature dependence through the Polyakov field is sufficient. 

The relevant order-like parameters will be analyzed and a special attention will be payed 
to the phase structure underlying the confinement- deconfinement and chiral transitions, as 
well as to thermodynamic quantities, such as the pressure, the behavior of the topological 
susceptibility and the convergence of meson axial and chiral partners. Section II is devoted 
to the description of the models, section III to the discussion of the results and, finally, we 
summarize our findings in section IV. 



II. MODEL AND FORMALISM 

We will use the SU(2) PNJL model with a 't Hooft interaction simulating the Ua(1) 
anomaly (instanton effects). Afterwards, several mechanisms of temperature dependence of 
the coupling coefficients are included, in particular the dependence of temperature through 
the Polyakov loop field, leading to entanglement of deconfinement and restoration of chiral 
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symmetry (EPNJL model). 

The PNJL Lagrangian with explicit chiral symmetry breaking where the quarks couple 
to a (spatially constant) temporal background gauge field (represented in term of Polyakov 



loops) is given by [3|, |23| : 



C vNjc = q(iY l D fl - m)q + £ 1 + £ 2 - U ($[A],$[A];T) , (1) 
with two different interacting parts 

A = 9i (W) 2 + (qil5Tq) 2 + (qfq) 2 + (qi^qf , (2) 

£2 = 92 {qqf + (qilsrq) 2 - {qfq) 2 - (qi-f 5 q) 2 . (3) 

The quark fields q = (u, d) are defined in Dirac and color fields, respectively with two flavors, 
Nf = 2 and three colors, N c = 3, the coupling coefficients g\ and g% have dimension energy" 2 , 
and rh = diag(m u ,m^) is the current quark mass matrix. For simplicity we assume below 
that m u = rrid = m. The original NJL model is reproduced for g\ = g 2 . 

The Lagrangian ([1]) is chiral invariant in the limit where the current quark masses vanish. 
Both terms C\ and £2 are invariant upon SU(2)i®SU(2)^®U(l) type transformations, but 
the £2 component makes the Lagrangian non-covariant upon Ua(1) transformations. The 



£2 term, that can be explicitly written in the form of a determinant (see [22|) and identified 
as an interaction induced by instantons, according to 't Hooft, explicitly breaks the axial 
symmetry even in the chiral limit. 

The quarks are coupled to the gauge sector via the covariant derivative = <9 M — iA^. 
The strong coupling constant gstrong has been absorbed in the definition of A 1 : A^{x) = 
9 strong A^(x)^f where A% is the SU C (3) gauge field and A a are the Gell-Mann matrices. 
Besides in the Polyakov gauge and at finite temperature A 11 = 5qA° = —i8%A . 
The Polyakov loop $ (the order parameter of Z 3 symmetric/broken phase transition in pure 
gauge) is the trace of the Polyakov line defined by: $ = -^-((Vexpi dr A 4 (x,t) )) „. 

The pure gauge sector is described by an effective potential U f$[A], $[A]; T) chosen to 
reproduce at the mean-field level the results obtained in lattice calculations: 

U ($ $-T) n(T) - 

{ ' ' ; = + 6(T)ln[l - 6$$ + 4($ 3 + $ 3 ) - 3($$) 2 ], (4) 



a cti a 2 63 



3.51 -2.47 15.2 -1.75 



TABLE I. Parameters for the effective potential in the pure gauge sector. 

where 



a (T) = a + ai (^j + a 2 and b(T) = b 3 f^j . (5) 

The effective potential exhibits the feature of a phase transition from color confinement 
(T < To, the minimum of the effective potential being at $ = 0) to color deconfinement 
(T > Tq, the minima of the effective potential occurring at $ ^ 0). 

The parameters of the effective potential U are given in Table [B These parameters have 
been fixed in order to reproduce the lattice data for the expectation value of the Polyakov 
loop and QCD thermodynamics in the pure gauge sector 24J, |25 |. 

The parameter T is the critical temperature for the deconfinement phase transition 
within a pure gauge approach: it was fixed to 270 MeV, according to lattice findings. This 
choice ensures an almost exact coincidence between chiral crossover and deconfinement at 
finite temperature, as observed in lattice calculations. 

The Lagrangian density [1] can be rewritten as 

Cpnjl = q(il fl - m)q + G s [{qqf + {qi^rq) 2 ] + G a [(qfq) 2 + (qi lb q) 2 ] 

-W($[A],$[A];T), (6) 

where, as quoted below, G s = g\ + g 2 is related to the n and a mesons and G a = g\ — g 2 to 
T) and ao mesons. 

The PNJL grand canonical potential density in the SU/(2) sector can be written as jsiH]: 



n^,^,M;T, l i)=U(^^,T)+AGN f (q l q l ) 2 -2N c N f / j^E, 

-2N f Tj ^3 (4(Ek) + zi(Ek)), (7) 

where Ej is the quasi-particle energy for the quark i: Ei — a/p 2 + M 2 , and and z J are 
the partition function densities. 
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The explicit expression of z$ and z$ are given by: 

z+(Ei) = Tr c ln [l + L^ Ei+ ^] = In {l + 3 ($ + $ e -^+")) e -/»(*+/0 + e -3/H^+M)} ; 

(8) 

= Tr c ln [1 + Le-^-ri] = In {l + 3 ($ + $ e -^-")) e -^-"> + e -W-f)} . 

(9) 

A word is in order to describe the role of the Polyakov loop in the present model. Almost 
all physical consequences of the coupling of quarks to the background gauge field stem from 
the fact that in the expression of 2$, $ or $ appear only as a factor of the one- or two- 
quarks (or antiquarks) Boltzmann factor, for example e~^^ Ei ~^ and e~ 2 ^^ Ei ~^' . Hence when 
$, $ — > (signaling what we designate as the "confined phase") only e -3 ^^* - ' 4 ' remains 
in the expression of the grand canonical potential, leading to a thermal bath with a small 
quark density. On the contrary, $,$—>• 1 (in the "deconfined phase") gives a thermal bath 
with all 1-, 2- and 3-particle contributions and a significant quark density H]- 

This formalism, presented here for completeness in the grand canonical approach, will be 
employed in the present work with /x = 0. This condition implies $ = $. 

We can redefine the coupling constants such as the set (g±, g 2 ) or (G s , G a ) will be replaced 
by (G, c) in the following parametrization: 

G 8 = gi + g 2 = G, G a = 9l -g 2 = G(l - 2c), (10) 



where c G {0, 1} 28|, |29j is a parameter that now specifies the degree of 11^(1) symmetry 



breaking. In the present work we take c = 0.2, a value which gives an adequate degree of 



anomaly in the vacuum [22[ , as will be shown in the sequel by the choice of model parameters. 
(Notice that g\ = G (1 — c) is the coupling constant of the four quark vertex associated with 
chiral symmetry effects, while g 2 = cG is the anomaly coefficient and is, for the present 
choice, 25% of g±). 

As it is well known, in pure gauge theory, the Polyakov potential induces a first-order 
phase transition at T = To. The PNJL model with the original To = 270 MeV, that 
reproduces pure gauge lattice QCD data, yields to a small difference between chiral and 
deconfinement transition temperatures. This value is however significant when we rescale 



T to 210 MeV, derived by RG considerations [30j so as to reproduce the lattice QCD 



result, Trf = 177 MeV for the deconfinement transition temperature as will be seen in 



the next section. Consequently, the PNJL result is not consistent with lattice QCD data 
for the transition temperatures. This indicates that the entanglement between chiral and 
deconfmement transitions is weak in this model. Following the suggestions of other authors 
3, |8J , we implement an explicit dependence of G on $ assuming the following form: 

= G [l - ai - a 2 ($ 3 + $ 3 )] , (11) 

which respect chiral, P, C and the extended Z3 symmetries. We use the parameters of 
reference 7j, where they have been fixed to reproduce the available lattice QCD data 
This leads to the values a\ = a 2 = 0.2 and T = 170 MeV. 

A standard calculation leads straightforwardly to the gap equation and to the meson 
propagators. The following gap equations are obtained: 

M t = m t - AG (qq). , (12) 

where one identifies i = u,d and Mj as the constituent quark mass. The quark condensates 
are determined by 

= -^r^ry = -iTr^p), (13) 

being Si(p) = (p 1 — Mi + ie)~ l the propagator of quarks. 

The mass spectra of the mesons is obtained by the analysis of the pole structure of the 
meson propagator, given by 

l-4G s<a U M (q 2 = M 2 M ) = 0, (14) 

where 

IW) = i I 0- 4 Tr [0 M S(p + q)0 M S(p)} (15) 

is the polarization operator for the quark-antiquark system regarding the channel with quan- 
tum numbers {M.} in the mesonic sector. As mentioned above, G s is related to ir and a 
mesons and G a to 77 and ao mesons. 

The topological susceptibility, \i is an essential parameter for the study of the problem 
of breaking and restoration of the Ua(1) symmetry. The topological susceptibility is defined 
as: 

X = j d 4 x (0|TQ(:r)Q(0)|0> c , (16) 
where c means connected diagrams, T the time order operator and Q(x) is of the form: 
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u 

[MeV] 


(gg)V3 
[MeV] 


m w 
[MeV] 


m a 
[MeV] 


m v 
[MeV] 


m ao 
[MeV] 


x 1/4 

[MeV] 


Model 


93 


-241 


140.2 


803.7 


704.5 


919.8 


180.8 


Exp. /Latt. 


92.4 


-267 


135.0 


400-1200 


547.3 


984.7 


180 



TABLE II. Numerical values for the calculated observables compared with experimental and lattice 
results, obtained with A = 590 MeV, G S A 2 = GA 2 = 2.435, c = 0.2 and m = 6 MeV. 



Q(x) = 2g 2 det [(g(l - 75 )g - det [(g(l + 7s ) 



(17) 



Taking only into account the connected diagrams, and only the terms of order 1/N C , fol- 
lowing a similar approach to 18|, we arrive at the following expression for the topological 
susceptibility (see |22j for details): 



X 



(18) 



where I\ = — {qq)jAMi. 

The present PNJL model has four parameters in the NJL sector: m, A, g\ and g 2 . We 
choose to adjust the parameters in vacuum by fitting to well known experimental data 
or lattice values: the mass of the pion, its decay constant, the quark condensate and the 
topological susceptibility that are shown in Table [TTJ The masses of the a, rj and ao mesons 
come as outputs. This set of parameters is crucial to get the correct description of isentropic 



trajectories in the T — Y limit 32] 



III. RESULTS AND DISCUSSION 



In a previous work 



22| , from the analysis of the behavior of the the topological suscepti- 



bility and of meson axial partners, we concluded that effective restoration of axial symmetry 
could only be achieved when g 2 was taken as a decreasing function of temperature. By 
using g 2 {T) = (72(0) / (1 + exp((T — T )/10)) we got the full restoration of axial symmetry 
even with a finite cutoff, A, at finite temperature; this ansatz is interpreted as an explicit 
mechanism of instanton suppression. In the present work, we will enlarge the investigation 
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of possible temperature dependence of the coupling coefficients and its consequences for 
several observables. We consider two scenarios, that in the PNJL model are (see Table HTT1) : 



Scenario A — We will keep g\ and g 2 as independent parameters. At finite temperature 
we may allow g 2 to have an explicit dependence on temperature but g\ is kept constant. 



• Scenario B — We will use the redefinition of Eq. fflO|) allowing for 

9l = G(l-c), g 2 = Gc. (19) 

In the last scenario g\ and g 2 are not independent, but G s = G will be kept always constant; 
on the contrary, G a varies since 

c(T) = 0.2/(T), where f(T) = 1/(1 + exp ((T - T )/10)) . (20) 

In this case only the topological susceptibility and the r) and ao meson masses will be 
affected. The other quantities have the same behavior as when gi and g 2 are kept constant. 

In the EPNJL we have equivalent scenarios with the replacement: g± — > <7i(<&) and 
92 9 2 {&)- In both cases we allow T (the critical temperature for the deconfinement 
phase transition within a pure gauge approach, usually fixed to 270 MeV) to have several 
values and we discuss this effect. We will compare results obtained in the framework of the 
PNJL and EPNJL models for the characteristic temperatures, the pressure, the topological 
susceptibility and the masses of chiral and axial meson partners. Along the work, we will 



2l|, for 



always consider the cutoff A — > oo at finite temperature [21|, [27j . As explained in 
T > T e ff (T > T e ff is the temperature at which the symmetry dynamically broken is fully 
restored, that is Mj = m,), the condensates are set to zero. The use of the infinite cutoff 
allows the presence of high momentum quarks, ensuring that the pressure goes to the Stefan- 
Boltzmann limit and, as shown in 32[ gives a better description of several thermodynamic 
quantities. 

As it can be seen in Table IHH in the PNJL model T x and Td never coincide but are 
closer for higher values of T , so it is adequate to T ~ 270 MeV, the value used in pure 
gauge approach. We also present results for T = 210 MeV, the value derived by RG 
considerations. The transition temperatures are defined by the peaks in the susceptibilities 
of the chiral condensate, for the restoration of chiral symmetry, and of the Polyakov loop, 
for de confinement- deconfinement transition. 
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I 


II 

(Scenario A) 


III 

(Scenario B) 


PNJL 


9i, 92 


9i, 92{T) 


gi(T),g 2 (T) 


EPNJL 


p 2 ($) 


0,($,T) 


ffl ($,T), 52 ($,T) 



TABLE III. Scenarios of the temperature dependence of g\ and 52 in PNJL and EPNJL models: 
in column I there is no explicit dependence on the temperature, columns II and III correspond to 
the scenarios A and B. 



Scenario B 


To 




T d 


A 


T eff 




[MeV] 


[MeV] 


[MeV] 




[MeV] 


PNJL 


210 


215 


177 


18% 


~ 250 




270 


237 


219 


8% 


~ 300 


EPNJL 


170 


173 


173 




~ 200 




270 


223 


223 




~ 300 



TABLE IV. Characteristic temperatures in PNJL and EPNJL model for different values of To 
(A = (T X -T d )/T x ). 



However, in the EPNJL model, where A = (T x — T d )/T x = by construction, we have 
more freedom to fix To and a lower value of To is convenient (~ 170 MeV) since it allows 
to reproduce lattice results for the critical temperature for deconfinement and restoration 



of chiral symmetry jfj]. The results are presented for scenario B, for reasons that will be 
explained latter. 

In the EPNJL model the coupling constants are replaced by effective couplings dependent 
on $, which ensures the entanglement between deconfinement and restoration of chiral 
symmetry. Now, using a EPNJL model with anomaly, a question arises: is the restoration 
of axial symmetry also entangled with two phase transitions mentioned above or should we 
have, like in PNJL model, an independent mechanism of instanton suppression? Having in 



mind the PNJL results [22] , we will discuss results for the following cases that are summarized 
in Table ED 

• G($), with <7i($), fi , 2( < ^ ) ), without extra dependence on temperature. 
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FIG. 1. Comparison of several quantities in the PNJL model (dashed lines) and EPNJL model (full 
lines), without explicit dependence of g\ and g 2 on temperature and at Tq = 270 MeV: coupling 
vertices (a), pressure (b) and meson masses (c). 

• G($), with <7i($,T), g 2 (§,T) - scenario B, the explicit dependence on temperature 
being introduced, as in PNJL model, through Eq. (E0|) . 

Lets us concentrate on observables related with the restoration of axial symmetry. Con- 
cerning the topological susceptibility, its vanishing is guaranteed in both models since the 
infinite cutoff leads to the vanishing of the quark condensate and x is proportional to (qq) 2 
(see Eq. (j!3p ). As for the gap between the masses of the meson axial partners, in order to 
get its vanishing one should have G a — >■ G s . We begin by considering no explicit mechanism 
of instanton suppression in the models, that is: we have g\ and g 2 constants, in the PNJL, 
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FIG. 2. Pressure for different values of To with scenario A in PNJL and EPNJL models (panels 
(a) and (b)); with scenario B in PNJL and EPNJL models (panel (c)). 



and gi(5>) and #2 (•=!?) in EPNJL (column I of Table III). As it can be seen in Fig. [U (a) for 
T = 270 MeV, although in the EPNJL model the effective vertices have a temperature de- 
pendence through the Polyakov field, $, G s does not converge to G a and, consequently, the 
axial partners (vr, a ) and (a, rj) do not degenerate in Fig. (JTJ (c)). Therefore, we conclude 
that in both models we need an explicit temperature dependence of the coupling coefficients 
in order to have effective restoration of axial symmetry. Finally, we notice that the results 
for the pressure are qualitatively similar in both models (Fig. HJ (b)). 

In order to check whether the scenarios explored for the temperature dependence are 
physically meaningful, we begin by plotting in Fig. [2] the pressure for different scenarios. 
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We see that in scenario A, for both models, there is a unphysical region of negative pressures 
for low values of To ((a) and (b)). This is probably due to a too fast instanton suppression 
when T takes lower values, so we need a mechanism to balance this effect if we want to 
consider low values of T . This problem does not exist in scenario B (Fig. [2J (c)), and, as will 
be seen (Figs. [3]and|I|), this ansatz also guarantees the restoration of axial symmetry. So, 
scenario B, which allows freedom in fixing To, will be the scenario adopted in the remaining 
of the present work. The difference in the behavior of the pressure in both scenarios can 
be understood from Eq. (J7|), which depends explicitly on the coupling constants through 
G s = gi + <72- The coupling G s that depends on c(T) (Eq. (ED}) in scenario A, but not in 
scenario B. 

In Fig. E]we plot the PNJL results for the meson masses, considering two value of T , 
in scenario B, and the topological susceptibility for different cases. We see that both the 
topological susceptibility as well as the gap between chiral and axial partners vanish, so axial 
symmetry is effectively restored. For T = 270 MeV the convergence of chiral partners occurs 
first than that of axial partners, as usual, with T e ff rj 300 MeV (T e ff is the temperature 
at which the effective restoration of the two symmetries is achieved). The new finding is 
that for To = 210 MeV the chiral and axial partners get degenerate very closely and T e ff 
is lower. An exploration of this temperature region could lead to a discussion about the 
sequence of restoration of two symmetries or even of its possible coincidence. Concerning 
the topological susceptibility fl3j (c)), the behavior is quite similar for scenario A and B 
and it vanishes at T « 325 MeV; when both g\ and g<i are kept constant, the topological 
susceptibility vanishes later, in this case this effect being only due to the vanishing of the 
quark condensate, a consequence of the full restoration of chiral symmetry. 

The behavior of the meson masses and topological susceptibility in the EPNJL model is 
qualitatively similar to PNJL model, as it can be seen in Fig. HI Here, also the restoration 
of chiral and axial symmetries become closer for low values of To, but the temperatures for 
the effective restoration of symmetries are slightly lower. We conclude that in EPNJL there 
is entanglement between restoration of chiral symmetry and deconfinement but not with 
restoration of axial symmetry. In order to have restoration of this symmetry an explicit 
correlated dependence on T of g\(§) and gi(<3>),that ensures instanton suppression and gives 
meaningful physical results, is necessary. 
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T (MeV) 

FIG. 3. Meson masses in the PNJL model with scenario B for two values of To and topological 
susceptibility with different scenarios. 

IV. CONCLUSIONS 

We performed a comparative study of deconfinement, restoration of chiral and axial sym- 
metries within the framework of the PNJL and EPNJL models. In the EPNJL model the 
coupling coefficients are endowed with a dependence on the Polyakov field, $, which guaran- 
tees that deconfinement and restoration of chiral symmetry occurs at the same temperature, 
which can be fitted to lattice result by choosing a low value for the parameter T = 170 MeV. 
In order to discuss the effective restoration of chiral and axial symmetries, we calculate the 
topological susceptibility and the masses of the mesons which are chiral and axial partners. 
At finite temperature we consider the cutoff A — > oo, in order to get the Stefan-Boltzmann 
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FIG. 4. Meson masses in the EPNJL model with scenario B for two values of To (panels (a) and 
(b)) and topological susceptibility with two scenarios (panel (c)). 

limit for the pressure. 

In both models we verified that the masses of chiral partners a and 7r converge and 
that the topological susceptibility vanishes, as a consequence of the effective restoration of 
chiral symmetry that occurs when the quark condensate vanishes. However, if no explicit 
temperature dependence of the coupling coefficient g 2 is considered, the masses of the axial 
partners do not converge and anomaly effects remain in the chirally symmetric vacuum. In 
order to get the vanishing of the anomaly it is enough that only g 2 is a decreasing function of 
the temperature that eventually goes to zero. However, if there is no mechanism to balance 
a fast decreasing of g 2 , for low values of T , an unphysical region of negative pressures 

1(3 



occurs. To prevent this, and have more freedom in the choice of To, we consider both g\ 
and (?2 as explicit functions of temperature, #2 going to zero while g\ attains a maximum 
value. This approach also guarantees the convergence of the axial partners, (it, a ) and 
{<J,f]). In both models the restoration of both symmetries is closer when low values of T 
are used. Our results follow the tendency of recent lattice results (calculated in the chiral 
imit), in particular in what concerns the degeneracy of mesons chiral and axial partners 
9), [it], 12]. Finally, we remark that although in the EPNJL model there is entanglement 
between restoration of chiral symmetry and deconfmememt, the restoration of the axial 
symmetry needs and extra assumption that ensures the annihilation of the contribution of 
the 't Hooft interaction. It should be noticed that the critical temperature, T x , signals only 
the partial restoration of chiral symmetry, its effective restoration occurring only at T e ff 
when the chiral partners get degenerate. So, even when chiral and axial symmetry become 
effectively restored, since T e ff > T x , some effects of the breaking of both symmetries remain 
above the critical temperature, in the interval (T x ,T e ff). 
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